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Abstract:  Final  Report 


Grant  #  USAF-AFOSR-79-0121 
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-r*“  Research  on  this  grant  has  focused  on  plasma  focus  experiments  in 
the  areas  of  particle  beam  generation  and  as  a  potential  repetitive 
opening  switch. 

In  pursuing  the  former  unique  diagnostic  tools  were  developed  to/, 
measure  the  scaling  of  particle  beam  current  and  energy  for  both  the 
electron  and  ion  beams  generated  by  the  device.  Simultaneous  measure-^ 
ments  of  the  energy  spectra  for  doth  the  electrons  and  ion  beams  were 
measured  for  the  first  time  as  were  scaling  laws  for  the  increase  of  ^ 
electron  energy  and  current  with  input  energy.  V. 

The  potential  of  the  plasma  focus  as  an  opening  switch  was  then  J 

n 

investigated.  Measurements  of  the  current  and  voltage  waveforms  indi-^j 
cated  that  the  resistance  of  the  pinch  was  roughly  ten  times  the  class'' 
cal  value  which  was  estimated  from  electron  temperature  measurements  and 
streak  pictures.  -Logd  experiments  were  performed  with  100MW  being  deJv 


vered  to  a  0.77  ohm  load 


a  pulse  compression  factor  of  10  and 


transfer  efficiency  of  2*10  .To  increase  the  efficiency  the  Impaler'j 
-  concept  was  devised  which  could  have  a  transfer  effi  of  well  over  50Z  , 
according  to  the  results  of  a  physical  model.  The  frequency  of  the  | 
microwave  emission  was  measured  using  the  delay  line  technique.  The  v 
observed  frequencies  were  most  consistent  with  the  lower  hybrid  fre-  T 

*  A, 

quency. ^A  model  is  being  devised  to  explain  this  emission  and  the  j 
enhanced  plfcnch  resistance. 
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Final  Report  for  AFOSR 
Grant#  USAF-AF0SR-79-0121 

1979-82  Restrike  Beam  Experiments  on  a  Dense  Plasma  Focus 

1982-84  Opening  Switch  Research  on  a  Plasma  Focus 


T  30  Sent.  1979  to  May  1982:  Restrike  Beam  Experiment 


The  goal  of  the  research  at  the  beginning  of  the  funding  period  was  to 
determine  the  nature  of  the  particle  beams  generated  by  a  dense  plasma  focus. 

This  research  was  performed  between  30  Sept.  1979  and  May  1982  when  results  of 
the  beam-target  model  for  neutron  production  generated  by  deuteron  beams  were 
presented  at  the  IEEE  Conference  in  Ottawa,  Canada.1  In  the  course  of  this 
research  unique  diagnostic  tools  were  developed  to  study  these  particle  beams 
Including  a  differentially  pumped  Faraday  cup  to  measure  ion  energies  by  time  of 

flight,  a  solid  state  nuclear  track  detector  technique  to  measure  high  fluences 

3  4  2 

of  light  ions,  an  electron  magnetic  spectrometer,  and  filtered  electron 

Faraday  cup.  These  and  other  diagnostic  tools  were  used  to  measure  the  energy 

spectra  and  intensity  of  the  electron  and  ion  beams  generated  by  a  plasma  focus 

and  the  scaling  of  these  beams  with  input  energy.  The  results  were  the  first  of 

their  kind  and  are  well  summarized  in  detail  in  our  publication  in  Nuclear 

Fusion5  which  is  included  here  as  an  appendix. 

While  the  scaling  laws  Indicated  a  hardening  (increased  particle  energy) 
and  intensifying  of  the  particle  beams  with  input  energy,  the  broad  spread  of 
particle  energies  (power  law  distributions)  seemed  to  limit  the  application  of 
such  an  accelerator  (although  the  electron  beam  may  have  promise  in  microwave 
generation5  or  as  a  pump  for  lasing  media.^).  At  the  same  time  it  became 
apparent  that  the  strong  current  interruption  in  the  plasma  focus  might  make  it  a 
good  candidate  as  an  opening  switch. 


.r.f  «r; 
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good  candidate  as  an  opening  switch. 

II.  May  1982  to  20  Sept.  1982 

This  research  began  after  meeting  Prof.  Schoenbach  at  the  3rd  Pulsed  Power 
Conference  in  Albuquerque,  NM  in  June  1981.  After  having  had  several  telephone 
conversations  with  him  since  that  meeting  it  was  realized  that  the  possibility 
of  a  sucessful  collaboration  existed.  This  was  strongly  encouraged  by  the 
grant  monitor  at  that  time,  Lt.  Col.  A.  K.  Hyder,  so  in  January  1982,  G.  Gerdin 
(principle  investigator  on  this  grant)  traveled  to  Texas  Tech  University  for  a 
three  day  visit  to  discuss  the  possible  nature  of  this  collaboration. 

Prof.  Schoenbach  is  an  ideal  man  with  whom  to  collaborate  since  he  had  been 
involved  with  dense  plasma  focus  (DPF)  research  for  about  ten  years  and  is  the 
author  of  over  ten  papers  on  the  subject.  Presently  he  is  directing  a  project 
on  diffuse  dischage  opening  switches  and  so  has  experience  in  both  DPF  and 
opening  switch  physics. 

A  previous  unsuccessful  attempt  to  use  the  plasma  focus  as  an  opening  switch 

was  regarded  as  not  a  true  test  of  a  plasma  focus  as  an  opening  switch.  That  is, 

it  did  not  appear  the  device  had  been  operated  in  an  optimum  manner  since  no 

o 

capacitive  discharge  data  showing  strong  interruptions  were  presented  and  the 

*  9  9 

initial  Iq  was  very  low.  Furthermore  the  circuit  used  represented  a  brute 

force  approach  and  the  results  indicated  a  more  sophisticated  circuit  or  gas 
handling  technique  would  be  required.  So  a  research  program  was  devised  to  pro¬ 
vide  experimental  tests  to  gain  an  understanding  of  the  DPF  interruption  physics 
under  ideal  conditions  of  a  capacitive  energy  drive  and  static  gas  fill. 

The  research  program  devised  on  this  trip  was  to  answer  these  four  basic 
questions: 


a)  Is  the  interruption  inductive  (L)  or  resistive  in  nature? 


b)  Can  the  power  be  tapped  into  a  load? 

c)  Can  the  timing  of  the  interruption  be  controlled? 

c)  Can  the  operation  be  made  repetitive. 

To  answer  the  first  question  it  is  important  to  know  the  nature  of  the 
resistivity  of  the  plasma  (classical  or  anomalous?).  To  determine  whether  the 
plasma  exhibits  classical  resistivity10  or  not  one  must  know  the  electron  tem¬ 
perature.  A  soft-x-ray  spectrometer  was  developed11  during  May-October  1982 

12—15 

which  utilized  pin  diodes  and  Ross  filters  for  fast  time  response.  Also  a 

computer  program  was  developed11  using  a  steady  state  corona  model. 

The  results  show  good  consistency,  yielding  temperatures  of  about  0.3  to 
0.5  keV,  regardless  of  the  type  of  impurity  gas  present  (neon  or  nitrogen; 

Figure  1).  This  result  is  encouraging  since  most  of  the  detected  x-ray  emission 
from  a  nitrogen  seeded  plasma,  is  of  the  continum  type,  as  opposed  to  a  neon 
seeded  plasma  where  line  radiation  dominates  the  detected  emission.  Roughly  the 
same  temperatures  are  expected  in  both  cases,  given  the  small  amount  of  impurity 
gas  present  (about  2%).  Strong  pinches  were  observed  with  pure  neon  as  the 
working  gas  (0.6  torr)  and  with  temperatures  about  the  same  as  reported  above. 

It  was  also  possible  to  calculate  the  classical  (Spitzer)  resistance  of  the 
pinch  from  evaluated  dimension  of  the  x-ray  image  of  the  pinch.  This  is  roughly 
9  nl)  at  a  temperature  of  500  eV.  This  value  is  about  fifty  times  lower  than  the 
observed  0.50  value  (from  current  interruption  measurements)  suggesting  the 
possibility  of  non-classical  (anomalous)  resistivity  playing  an  important  role 
in  the  focus  current  interruption  at  pinch  time. 
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Figure  1.  Ratio  of  band  strengths  vs,  temperature,  as  determined  by  computer 
simulation  (corona  model)  for  deuterium  fills  with  2%  nitrogen  and 
2%  neon. 
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During  this  period  the  research  effort  on  the  plasma  focus  as  an  opening 
switch  benefited  greatly  by  the  collaboration  with  the  pulsed  power  group  at  the 
Texas  Tech  University.  This  collaboration  took  the  form  of  a  winter  visit  by 
Professors  Schoenbach  and  Krompholz  to  the  University  of  Illinois  to  plan  the 
preparations  for  Professor  Krompholz's  more  extended  visit  in  June.  During  his 
visit  experiments  involving  microwave  emission,  pinch  phenomena  through  streak 
pictures  and  the  Mark  I  load  experiments  were  performed. 

The  research  could  be  divided  into  three  phases: 

1)  Development  and  testing  of  a  0«D  model  describing  the  motion  of  the 
current  sheaths  in  a  dense  plasma  focus. 

2)  Observation  of  the  phenomena  occurring  at  pinch  time  with  normal 
operating  procedures . 

3)  Performing  DPF  load-coupling  experiments. 

2  3 

The  0*D  model  developed  was  a  combination  of  a  couple  of  others  *  and  was 

found  to  be  in  good  agreement  with  observations  of  the  motion  of  the  current 

u 

sheath  in  the  pre-pinch  stage.  The  model  predicted  that  the  pinch  would  bounce 
radially  and  result  in  a  rapid  voltage  reversal  across  the  load  and  a  loss  of 
power  flow.  This  prediction  was  tested  by  the  experiments  of  phase  2  and  phase 
3. 

Phase  2  of  the  experiments  was  the  major  activity  during  the  first  two 
weeks  of  Prof.  Krompholz's  visit.  These  experiments  consited  of  the  following: 

1)  Observation  of  the  time  behavior  of  the  microwave  emission  over  various 
frequency  bands. 

2)  The  use  of  radial  and  axial  streak  camera  pictures  to  observe  phenomena 
at  pinch  time  to  test  some  of  the  conclusions  from  phase  1. 


3)  Taking  soft-xray  pinhole  camera  pictures  simultaneously  with  the  axial 
streak  picture  to  get  a  more  comprehensive  view  of  the  pinch  phenomena. 

Microwave  emission  was  observed  simply  by  placing  microwave  horns  and 
detectors  in  the  room  outside  the  device.  Emission  was  observed  to  be  coming 
from  between  the  parallel  plates.  The  period  of  duration  of  the  emission  gets 
shorter  as  the  frequency  of  the  band  is  increased  (Figure  2)  with  the  highest 
frequency  emission  occuring  just  after  pinch  time.  This  could  be  consistent 
with  an  emission  frequency  proportional  to  the  local  magnetic  field  (e.g.  the 
lower  hybrid)  frequency  or  density  such  as  the  ion  plasma  frequency.  The  magni¬ 
tude  of  the  emission  frequency  is  most  consistent  with  the  lower  hybrid  for  con¬ 
ditions  believed  to  be  present  in  the  discharge  during  this  period  (see  next 
section).  Since  the  lower  hybrid  instability  is  associated  with  anomalous 
resistivity  this  could  prove  to  be  an  extremely  interesting  result. 

The  radial  streak  pictures  provided  a  test  as  to  whether  the  radial  bounce 

predicted  by  the  0*D  code  actually  occurred.  The  experiments  with  a  streak 

camera  set  up  to  observe  the  radial  implosion  of  the  pinch  (Figure  3)  confirmed 

the  presence  of  this  bounce;  the  effect  of  the  bounce  limiting  the  power  flow  to 

1 9 

the  load  was  confirmed  by  the  results  of  the  load  coupling  experiments. 

Axial  streak  pictures  (Figures  4,  5)  revealed  that  the  first  compression 
involves  an  approximately  uniform  compression  of  the  plasma  column  but  that 
second  compression  involves  local  constrictions  of  the  column  indicating  the 
onset  of  mhd  instabilities  (usually  kinks).  Comparison  with  the  current  wave¬ 
form  of  the  DPF  (Figure  5)  indicates  that  the  maximum  device  current  interrup¬ 
tion  occurs  during  1st  compression. 

Thus,  the  conclusions  resulting  from  phase  2  of  the  experiments  were  that 
the  radial  bounce  occurs  and  is  detrimental  to  efficiency  energy  transfer  and 
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Figure  2.  Typical  microwave  power  generated  by  the  Illinois  plasma  focus 
(not  the  same  shot)  as  observed  with  a  horn-detector  and 
oscilloscope  over  the  following  frequency  bands:  1)  G  band 
(4-6  GHz),  2)  X  band  (8-12  GHz),  3)  P  band  (12.4-18  GHz),  and 
4)  R  band  (26.5-40  GHz). 
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Figure  3.  Top:  A  schematic  of  the  plasma  and  streak  camera  in 
the  radial  streak  configuration. 

Middle:  Camera  view  of  the  plasma  at  three  successive 
instants  of  time. 

Bottom:  An  interpretation  of  the  radial  streak  picture 
and  an  actual  photograph. 
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Axial  Streak  Pictures 
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Figure  4.  Top:  A  schematic  of  the  plasma  and  streak  camera  in  the  axial 
configuration.  Middle:  Various  possible  macroscopic  or  mhd 
Instabilities  that  could  occur.  Bottom:  The  appearance  on 
the  screen  if  either  type  instability  (sausage  or  kink)  were  to 
occur. 
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Current  Waveforms  and  Axial  Streak  Pictures 
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Figure  5.  Comparison  of  an  axial  streak  picture  with  the 
time  derivative  of  the  device  current.  The  Deak 
in  the  current  interruption  occurs  during  first 
compression. 
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that  first  compression  Is  the  most  Important  phase  of  energy  transfer.  Thus, 

reducing  the  effects  of  the  Instabilities  and  enhancing  the  effects  of  the  first 

compression  appear  to  be  the  keys  to  Improving  efficiency. 

The  DPF  load  coupling  experiments  (phase  3)  tended  to  confirm  the  results 

of  phase  2.  The  load  was  hooked  up  to  the  DPF  across  the  breach  end  of  the 

device.  The  results  are  shown  In  Table.  I.  Coupling  efficiency  appears 

limited  by  rapid  oscillations  In  load  voltage  (Figure  6)  which  appear  to  be 

caused  by  the  radial  bounce  and  subsequent  mhd  activity  during  2nd  compression 

(Figure  5).  The  voltage  across  the  load  at  the  breach  end  of  the  DPF  Is  con- 

1 8 

sis tent  with  the  0-D  computer  model. 

One  strategy  to  Improve  the  coupling  efficiency  including  the  addition  of 

an  axial  magnetic  field  to  stabilize  the  column  to  enhance  first  compression  and 

12 

prevent  mhd  instabilities.  The  addition  of  an  axial  magnetic  filed  has  been 
1 2 

previously  found  to  affect  the  discharge  as  follows: 

1)  Stabilize  the  plasma  column  to  mhd  instabilities 

2)  Increase  the  duration  of  current  interruption 

3)  Decrease  the  magnitude  of  the  rate  of  current  Interruption. 

20 

The  first  result  is  consistent  with  classical  mhd  theory  and  the  last  two  have 

19 

been  confirmed  by  the  results  of  preliminary  experiments  on  the  Illinois  DPF. 

The  third  result  Implies  a  reduced  load  voltage  which  hurts  coupling  efficiency 
and  thus  must  be  considered  as  a  trade-off  in  the  design  of  an  optimal  energy 
coupling  device. 

In  conclusion  the  results  of  the  FY  1982-83  research  on  this  project  have 
supplied  at  least  partial  answers  to  the  first  two  questions  listed  in  the 
research  objectives  sections  of  this  report.  That  is  the  current  interruption 
is  both  resistive  and  inductive  in  nature  as  evidenced  by  partial  reversal  of 


VV 


t 


12 


TABLE  1 


Mark  II 

Impale r 

Bank  Capacitance,  C 

40pF 

40pF 

Bank  Voltage,  V 

o 

25kV 

25kV 

Bank  Energy,  E 

o 

12.5kJ 

12.5kJ 

Plasma  Current,  1 

560kA 

560kA 

Load  Inductance,  L^ 

30nH 

26nH 

Load  Resistance, 

0.770 

0.77£1 

Time  to  Interruption, 

l*75|is 

1 .75ps 

Period  of  Interruption,  T 

200ns 

100ns 

Peak  Load  Current,  1^ 

1  lkA 

250kA 

Peak  Power  to  Load,  P^ 

100MW 

63GW 

Peak  Load  Voltage,  V^ 

17.5kV 

250kV 

Energy  Transfer, 

25J 

6.3kJ 

Transfer  Efficiency,  n 

2*10-3 

50% 

Load  Experiment*. 
Rl  =  0.77  £1 
Ll=30  nH 
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Figure  6.  Experimental  results  from  the  Mark  II  load 
coupling  experiment  (Table  I). with  a  con¬ 
figuration  similar  to  Figure  5  except  that 
load  switch  S2  was  closed  at  all  times;  hence 
the  occurrence  of  load  current  before  pinch 
time.  The  energy  transferred  during  the  first 
sharp  load  voltage  spike  was  27  J. 


load  voltage  due  to  the  radial  bounce  which  Is  an  Inductive  effect.  Also  power 


can  be  tapped  Into  a  load  but  enhancing  the  efficiency  to  an  Interesting  power 
transfer  device  will  probably  involve  strategies  to  increase  the  resistive 
nature  and  suppress  the  inductive  nature  (IL)  of  the  device. 

IV.  30  Sept.  1983  to  30  Sept.  1984 

A)  Research  Objectives 

In  the  previous  year,  100  MW  was  coupled  into  a  load  but  the  efficiency 
(energy  delivered  to  the  initial  capacitor  bank  energy)  was  only  0.2%.  Thus  the 
research  centered  on  discovering  ways  to  greatly  enhance  this  efficiency.  To 
achieve  this,  a  two  pronged  approach  was  formulated:  one,  to  try  to  understand 
the  nature  of  the  current  interruption  in  the  plasma  focus;  and  two,  to  use  cir¬ 
cuit  analysis  and  preliminary  measurements  to  design  a  potentially  more  effi¬ 
cient  configuration  for  testing.  The  progress  made  in  achieving  these  goals  is 
presented  in  the  next  section. 

B)  Status  of  the  Research  Effort 

The  research  performed  in  the  past  fiscal  year  involved  the  detailed  study 

of  microwave  emission  by  the  device  and  the  design  of  a  muzzle-end  load  coupling 

experiment  to  enhance  the  voltage  delivered  to  the  load.  Since  the  microwave 

2 1 

emission  is  orders  of  magnitude  above  thermal  levels,  it  is  an  indication  of 
22  23 

microturbulence  and  perhaps  enhanced  resistivity  .  If  the  latter  can  be 
monitored  and  controlled,  it  could  provide  the  basis  for  a  high  power  opening 
switch.  Muzzle-end  load  coupling  is  predicted  to  increase  the  load  voltage  by 
a  factor  of  five  greater  than  the  25kV  breach  voltage  pulse  we  observe.  Since 
the  power  to  resistive  and  hence  transmission  line  loads  is  proportional  to  the 
square  of  this  voltage,  the  system  efficiency  may  be  increased  by  a  factor  of 
twenty  five.  The  progress  in  the  muzzle-end  coupling  research  and  the  results 


This  work  consisted  of  two  phases:  the  first  involving  preliminary  experi¬ 
ments  and  the  second,  involving  the  coupling  apparatus,  was  designed  and 
constructed.  Preliminary  measurement  of  the  coupling  efficiency  will  be  made  in 
January  1985. 

To  get  an  idea  of  the  nature  of  the  potentials  in  the  muzzle  end  of  the 

device  a  floating  Langmuir  probe  was  inserted  along  the  two  lines  labeled  z  and 

r  in  Figure  7.  Since  the  impedance  of  the  floating  probe  is  kept  high,  the 

probe  acquired  a  potential  equivalent  to  that  of  a  muzzle-end  coupling  electrode 

before  it  starts  to  draw  significant  current.  Hence  the  floating  potential, 

25 

which  is  below  the  true  plasma  potential  by  3.87  •  KT  /e  for  a  deuterium 

e 

plasma  (~-100V  in  the  current  sheath  ~-l  kV  in  the  pinch),  is  really  the  quan¬ 
tity  of  interest  for  load  coupling. 

The  results  are  shown  in  Figure  8  and  typical  waveforms  in  Figure  9.  They 
are  somewhat  disappointing  in  that  in  no  case  is  the  floating  potential  greater 
than  the  breach  voltage  in  spite  of  the  predictions  of  the  model.  Either  the 
high  potential  region  is  very  localized  on  the  axis  and  missed  by  the  probe  (the 
probe  was  bent  when  placed  on  the  axis  and  hence  the  measurement  of  that  point  is 
unreliable)  or  the  model  is  not  valid. 

The  muzzel-end  coupling  device  was  then  designed  to  test  these  results  in  a 
practical  manner  and  to  test  a  novel  second  approach  to  current  interruption 
which  we  call  the  Impale r  concept . 

The  impaler  concept  is  to  achieve  a  current  interruption  by  'impaling'  the 
current  sheath  on  a  ridged  insulator  as  shown  in  Figure  10.  This  configuration 


Figure  8.  The  ratio  of  the  peak  in  the  floating  potential  (F.P.)  to  the 
peak  in  the  breach  voltage  (C.P.)  as  measured  along  the  r  and 
z  tracks  shown  in  Figure  7.  Along  the  r  tract  (top)  the  radial 
positions  of  the  outer  edge  of  the  center  electrode  (C.E.)  and 
the  inner  edge  of  the  outer  electrode  (O.E.)  are  shown.  Similarly 
the  position  of  the  tops  of  the  outer  and  center  electrodes  are 
indicated  along  the  z  track  (bottom). 
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between  the  center  and  load  electrodes  which  is  represented  by  the  switch  S9  in  Figure 


,VvV.  V.  VV  •  z  ■ 


utilizes  the  continuation  of  the  motion  of  the  current  sheath  from  the  rundown 
phase  to  drive  the  interruption.  The  plasma  itself  provides  the  load  switch 
connection  S2  in  the  opening  switch  circuit  (Figure  11)  and  by  adjusting  the 
relative  positions  of  the  load  electrode  and  the  impaler  insulator,  it  should  be 

possible  for  the  load  connection  and  the  current  interruption  to  be  made 

simultaneously.  Thus  prepulses  on  the  load  may  be  effectively  eliminated.  To 
study  muzzle-end  coupling  the  focus  without  the  impaler  effect  one  just  has  to 
place  the  impaler  insulator  further  away  from  the  plasma  focus  electrodes  with 
respect  to  the  load  electrode.  Then  the  pinch  occurs  before  the  impaling  of  the 
current  sheath  takes  place. 

The  scaling  laws  for  the  impaler  concept  are  quite  favorable  and  are  based 

1 2 

on  the  better  understood  rundown  phase  of  the  plasma  focus  discharge.  The 
interruption  time,  T,  should  be  equal  to  the  thickness  of  the  current  sheath, 

lg,  divided  by  its  speed  vg.  For  plasma  focus  devices  this  speed  is  given  by:16 


(  — £ -  ) 

2woPo 

o  o 


where  pQ  is  the  mass  density  of  the  filling  gas.  Thus  the  voltage  generated  by 


interruption,  V^,  is  : 


Vt  -  LgI  *  LgI  Vg / ^  (2) 

8 

where  L  and  I  are  the  gun  inductance  and  the  device  current  respectively.  For 

u 

1 8  •  '  7 

values  measured  in  our  experiment  (Z  ~  1cm,  v  ~  10  cm/sec,  L„  ~  50nH,  and 

8  8  U 

—7 

I~500kA),  T  and  V  are  10  sec  and  250  kV  respectively.  For  resistive  loads 
(and  hence  long  transmission  line  loads)  the  power  transfer,  P^,  is  V^/R^, 
where  is  load  resistance.  Therefore,  is  given  by: 


LrI  V8/,  2 


•>  „  ’,*•  ,•  "  V  •  V*  »’•  V*  /•  *.  *  *.  *•  ‘.N  *•  l’»  “•  *.  **,  *.  *•.  •, 
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storage  inductor  and  S,  is  the  opening 


where  p  is  the  original  pressure  of  the  fill  gas.  The  energy  to  the  load, 
is  then  given  by: 


E.  ~  P.T  * 
l  l 


T  2  T2 

lg  1 


V 


V 


h 
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The  system  efficiency,  n,  them  scales  as: 


n  “  Ei/. 


LGX  LG  Vs 


2LG  Vs 
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The  parameters  for  our  present  experiment  are  shown  in  Table  I  along  with 
our  previous  load  coupling  results. 

The  improvement  in  performance  is  incredible.  If  the  sheath  thickness 
doesn't  change  with  bank  current  one  might  be  able  to  achieve  a  transfer  effi¬ 
ciency  of  one  into  a  resistive  load  at  a  device  current  of  1.1  MA.  Such  currents 
are  relatively  common  in  other  plasma  focus  devices  and  we  could  even  achieve 
this  value  in  our  present  device  with  present  facilities  (although  the  spark 
gaps  would  have  to  be  configured). 

While  this  analysis  has  been  based  on  several  assumptions,  the  fact  that  the 


current's  own  magnetic  pressure  drives  the  current  interruption  makes  the  con¬ 
cept  promising.  The  impaler  electrodes  have  been  constructed  and  the  experi¬ 
ments  in  the  coming  months  could  be  very  interesting. 
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B.2  Microwave  Emission 

Due  to  the  possible  relationship  between  the  generation  of  microwave 
radiation  and  turbulence  induced  resistivity  of  the  plasma,  the  nature  of  the 
microwave  emission  was  studied  in  detail  to  test  various  models  for  these 
phenomena. 

The  emission  model  is  being  formulated  on  the  basis  of  various  emission 
experiments  which  generally  fall  into  three  categories.  They  are  as  follows: 

A) .  Measurement  of  source  location  and  source  strength  for  various 

microwave  bands 

B) .  Measurement  of  the  frequency  of  the  highest  (and  shortest  duration)  jf 

emitted  signal  by  the  time-of-f light  spectorscopy. 

C) .  Time  sequencing  of  the  emission  with  various  phenomenon  that 

characterize  a  plasma  focus  discharge. 


In  the  following  three  subsections  each  of  these  experiments  will  be 
discussed  as  to  the  nature  of  the  measurement  and  the  results.  On  the  basis  of 
the  results  of  experiments  (continuing)  model  will  be  selected  in  the  fourth 
subsection. 


B.2.1  Microwave  source  location  and  source  strength 


Microwave  emission  was  observed  by  placing  microwave  horns  and  detectors 


outside  the  device.  This  is  consistent  with  the  observation  from  the  DPF  devi- 
21  26 

ces  at  Darmstadt  and  Stevens  Institute.  The  power  levels  are  on  the  order 

2 

of  ~1  mW  for  the  R  band  (21  to  40  GHz)  (or  an  intensity  of  0.1  W/m  )  and  somewhat 


higher  for  P,  X,  and  G  bands  (Figure  2);  the  emission  occurs  during  pinch  time 

when  the  current  is  being  interrupted.  The  emission  is  polarized  in  the  direc— 

27 

tion  of  the  normal  of  the  parallel  plates  of  the  capcitor  bank. 


r 
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The  source  of  the  microwave  emission  may  be  determined  by  changing  the 
locations  of  the  horns  between  shots  and  noting  the  change  in  signal  strength. 
From  these  experiments,  the  microwaves  appear  to  be  emitted  from  between  the 
plates  of  the  parallel  plates  of  the  transmission  lines  and  capacitor  banks;  the 
strength  in  the  R  band  (21  to  40  GHz)  is  roughly  uniform  with  position.  This 
suggests  that  the  plates  are  acting  as  antennas  and  transmission  lines  as 
opposed  to  waveguides.  This  possibility  is  strongly  supported  by  the  obser¬ 
vation  of  emission  in  the  G  (4-6  GHz)  and  X  bands  (8-12  GHz)  during  periods  when 
shorter  wavelength  emission  in  the  P  band  (12-18  GHz)  and  the  R  hand  did  not 
occur.  The  cutoff  for  the  parallel  plates  of  the  capacitor  banks  (if  they  were 
to  act  like  waveguides)  is  14  GHz  and  no  emission  in  the  G  and  X  bands  should  be 
observed. 

To  obtain  more  information  on  the  location  of  the  microwave  source  an  R- 
band  detector  was  attached  to  the  anode  so  that  it  could  view  the  region  between 
the  inner  and  outer  electrode  directly  (Figure  12).  The  waveguide  was  aligned 
so  the  electric  field  vector  for  its  TE^  ^  mode  is  parallel  to  the  radial  vector 
of  the  coaxial  electrodes.  The  intensity  level  measured  was  typically  over  600 
W/m^  and  occured  during  the  interruption  in  current  (Figure  13).  This  intensity 
is  almost  104  higher  than  that  measured  outside  the  device.  Since  the  wave¬ 
guide  was  pointed  toward  the  muzzle  end  of  the  coaxial  electrodes  (Figure  12) 
this  is  strong  evidence  that  the  source  region  can  be  restricted  to  muzzle  end 
of  the  device. 

B.2.2  The  measurement  of  the  frequency  of  the  R-band  emission 

Microwave  radiation  in  the  R-band  is  typically  30  to  50  ns  in  duration  and 
contains  several  peaks  each  with  a  width  of  less  than  five  nanoseconds 
(FigurelA).  The  narrowness  of  the  peaks  makes  it  possible  to  measure  the  fre- 
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Figure  13.  Typical  waveforms  for  the  breach  emission  of  microwaves  in  the 

R  band.  Also  shown  is  a  current  waveform  which  shows  the  emission 
occurs  near  the  end  of  the  current  interruption  as  consistent  with 
previous  measurements.^ 
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quency  of  the  emission  by  the  delay  line  technique.  *  This  technique  utili¬ 
zes  the  fact  that  near  the  cut  off  frequency  of  a  waveguide,  f£,  the  group 
velocity,  v  of  an  electromagnetic  pulse  is  frequency  dependent  and  obeys  the 

o 

dispersion  relation 

vg  -  cg  0-  fc2/f2)'*  (7) 


where  c  and  f  are  the  speed  of  light  in  the  guide  and  the  frequency  of  the 

O 

pulse  respectively.  Thus  for  our  experimental  arrangement  shown  in  Figure  15, 
the  microwaves  traveling  down  the  7.35  meter  of  waveguide  will  experience  a  fre¬ 
quency  dependent  delay,  t,  with  respect  to  those  detected  by  the  10cm  length 
guide  for  the  other  detector.  The  delay  and  pulse  frequency  are  related  by 


f  -  f  (l-d2/cV)^ 


(8) 


where  d  and  c  are  the  difference  in  lenghts  of  the  two  guides  (7.25m)  and  the 

speed  of  light,  f  is  21.08  GH  for  the  TE,  n  mode. 

c  z  l  ,u 

Higher  order  modes  either  1)  will  not  propagate  in  the  parallel  plate 

transmission  lines  of  the  device  and  capacitor  bank  (TM),  2)  will  not  be  detected 

due  to  the  placement  of  the  antenna  in  the  system  detector  (TE  ;  m,n  even), 

m,n 

or  3)  have  a  cutoff  frequency  higher  than  the  maximum  frequency  to  which  the 
detector  will  respond  (TE^  q).  Since  the  polarity  of  the  emission  is  perpen¬ 
dicular  to  the  parallel  plates  (see  section  B.2.1)  and  the  TE.  .  has  a  higher 
attenuation  per  unit  length,  the  TE^  Q  is  the  mode  detected  by  this  method. 

The  relationship  between  the  delay  time  and  frequency  (Equation  8)  is  shown  in 
Figure  16. 

The  delay  line  technique  has  been  previously  used  to  measure  the  emission 


spectra  of  microwaves  generated  in  the  propagation  of  relativistic  electron 
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Figure  15.  Schematic  of  the  time  of  flight  spectrometer. 
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Figure  16.  Transit  time  of  a  microwave  pulse  through  7,22  m  of  R-band  (a=7,ll  mm,  b=3,56  mm)  waveguide 


beams  through  gases30  and  to  measure  the  output  frequency  of  an  electron 

29 

synchrotron  laser  driven  by  an  intense  relativistic  electron  beam.  In  addi- 

31  5 

tion,  it  was  used  to  study  the  emission  of  microwaves  by  the  electron  beam 
generated  by  a  plasma  focus.  These  microwaves  were  emitted  from  a  hole  in  the 
center  electrode  and  were  found  to  be  high-power,  narrowband  millimeter  waves 
(P~  1MW;  Af/f  <  5  •  10'2;  30  GHz  <  f  <  120  GHz). 

The  origin  of  these  latter  microwaves  is  not  the  same  as  that  of  the 
microwaves  studied  here  since  a  metal  vacuum  chamber  below  the  center  electrode 
prevented  observation  of  electron  beam  generated  microwaves  and  as  will  be 
demonstrated  in  the  next  section  the  period  of  emission  of  the  microwaves  does 
not  coinside  with  the  emission  of  the  electron  beams.  Furthermore  the  current 

i  C  3 

sheaths  (Figure  12)  are  of  sufficient  density  (>10  /cm3)  and  thickness  (~lcm) 

to  cutoff  any  millimeter  waves  from  propagating  from  the  region  of  the  hole  in 

the  center  electrode  to  the  region  between  the  inner  and  outer  electrodes 

(Figure  12);  no  microwave  radiation  was  detected  coming  from  anywhere  but  from 
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between  the  parallel  plate  of  the  capacitor  bank. 

The  experimental  arrangement  is  shown  in  Figure  13.  After  a  long  experi¬ 
mental  run  twelve  shots  were  found  to  have  sufficiently  few  peaks  in  the  R-band 
so  that  the  delayed  peak  and  its  corresponding  non-delayed  peak  could  be  iden¬ 
tified.  Two  examples  are  shown  in  Figure  17.  After  corrections  are  made  for 
varying  cable  lengths  the  frequency  of  the  pulse  may  be  determined. 

From  the  twelve  shots,  fifteen  frequencies  could  be  determined.  Two  fre¬ 
quencies  were  found  to  be  28  GHz  (±2GHz)  and  the  rest  lay  between  23  and  24  GHz 
(±1G1Iz).  Since  these  frequencies  are  close  to  cutoff  (Figure  13)  the  accuracy 
of  the  measurement  is  quite  good.  All  the  data  were  taken  with  a  capacitor  bank 
energy  of  12.5  kJ  and  a  deuterium  gas  pressure  between  3.0  and  3.6  torr. 


Since  it  has  been  hypothesized  that  this  emission  is  caused  by  current 
21 

driven  instabilities  it  is  of  interest  to  compare  these  measurements  with 

frequencies  of  the  dominant  modes  of  each  of  the  current  driven  instabili- 
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ties  which  may  be  candidates.  These  instabilities  are: 

1)  Buneman  (vD  >  v  h  ) 

2)  Ion  Acoustic  (T  »T .  ) 

e  l 

3)  Lower  hybrid  (B<1) 

4)  Anisotropy-driven  electromagnetic  (T^^t^j) 

where  vD,  v  ,  Te>  and  3  are  the  drift  velocity  of  the  electrons,  their 
thermal  velocity,  their  temperature,  that  of  the  ions,  and  the  ratio  of 
kinetic  to  magnetic  pressure  respectively.  The  plasma  parameters  for  the 
plasma  focus  during  the  period  of  microwave  emission  are  as  follows: 

Table  II 

Device  Current,  I  390kA 

Plasma  radius,  r  0.39cm 

Magnetic  field  at 

the  plasma  surface  B  20  Tesla 

Plasma  Temperature11*32  300  eV 

Plasma  density33  5«1018/cm  3 

The  corresponding  peak  frequencies  of  the  various  current  driven  instabili¬ 
ties  are  shown  in  the  following  table: 


-  *  *  w  •  V-  V*  V*  *  m  •«  *  j,  •  V  *  w  Vw 


Table  III 


Instability 
Buneman 
Ion  Acoustic 
Lower  Hybrid 


Maximum  Growth,, 
2  3 

Rate  Frequency 


l/2(me/2mi)l/3  f 


pe 


fpi  "  3 

9/4  fLR  =  9/4  (fcefci) 


V2 


Anisotropy-driven 
electromagnetic  f  ^ 


Estimate  of 
the  Magnitude 

520  GHz 

192  GHz 

21  GHz 


0.15  GHz 


where  fpg,  fpi,  fce>  ^ci*  an<*  ^LH  are  e^ectron  anc^  ion  plas™a»  electron  and 
ion-cyclotron  and  lower  hybrid  frequencies  respectively. 

The  Buneman  and  anisotropy-driven  electromagnetic  instabilities  have 
frequencies  which  differ  the  most  from  the  measured  value  and  hence  probably 

can  be  eliminated  from  consideration  as  source  candidates.  The  ion  acoustic 

16  3 

instability  could  act  as  a  source  in  low  density  regions  (n  <7*10  /cm  )  but 

6 

the  necessary  condition  for  it  to  exist,  Tg»T1  is  very  unlikely  to  occur  in 
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a  shock  heated  plasma. 

Hence  the  most  likely  candidate  is  the  lower  hybrid  instability  which 
occurs  in  regions  of  strong  magnetic  field  gradients  and  hence  the  value  quoted 
in  the  table  is  really  an  upper  limit  to  its  magnitude. 

B.2.3  Time  sequence  of  the  microwave  emission  with  various  phenomena  that 
characterize  a  plasma  focus  discharge. 


In  order  to  test  possible  physical  models  the  time  of  the  microwave 
emission  was  compared  with  several  other  phenomenon  associated  with  the 
discharge.  These  other  phenomenon  were: 


Phenomenon 


Diagnostic 


1)  Time  derivative  of  the  device  current 


Rogowski  coil 


2)  Soft  Xray  emission  Pindiode  detectors 

and  various  filters 

3)  Current  sheath  dynamics  Radial  and  axial 

streak  pictures 

4)  Hard  Xray  and  neutron  emission  Fast  scintillator 

and  PM  tube 

5)  Electron  beam  emission  in  various  energy  bands  Magnetic  Spectoroeter 

The  magnetic  spectrometer  was  attached  to  the  bottom  of  the  anode  so  it 
could  sample  the  electron  beam  passing  down  the  hole  in  the  center  electrode 
as  shown  in  Figure  18.  The  ranges  of  energy  sampled  by  a  given  channel  of  the 
spectrometer  are  given  in  Figure  19.  Also  shown  in  Figure  18  is  the  location  of 
the  Rogowski  coil  to  which  was  used  to  measure  the  time  derivative  of  the  device 
(or  main  bank)  current.  The  arrangement  to  measure  the  soft  Xray  emission  is 
shown  in  Figure  20;  the  long  collimating  tube  betweeen  the  pindiodes  and  the  DPF 
chamber  restricts  the  observed  region  of  the  pinch  to  a  spot  on  the  device  axis 
2ram  in  diameter  and  5mm  above  the  center  electrode.  The  photomultiplier  tube 
and  scintillator  were  placed  3.96m  away  from  the  position  of  the  pinch  and  the 
tube  has  inherent  delay  of  40  ±7ns.  Corrections  for  time-of-f light  and  cable 
delays  have  been  made  in  the  data  reported  below.  No  external  magnetic  field 
was  applied  to  the  discharge  during  the  experiments  reported  here. 

The  results  for  the  various  phenomena  excluding  the  streak  picture  are 
shown  in  Figure  21  and  typical  comparisons  between  signals  are  shown  in 
Figures  22-24.  The  two  soft  Xray  pulses  are  due  to  two  densif ications  along  the 
device  axis:  the  first  being  the  initial  compression  of  the  plasma  column  and 
the  second  several  seperate  densif ications  along  the  axis.  The  latter  signals 
the  breakup  of  the  pinch  and  is  probably  caused  by  mhd  instabilities(Figure  5). 

These  results  show  that  the  electron  beam  and  the  microwaves  are  not  directly 
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Figure  18.  Electron  Beam  Spectrometer. 


EVENT 


o 

<3-  LU 


4J 

XT*— 

U  S  ^ 
^  fOXH-* 
O  CJ3TJ 

u_ 

0  4-^0 
E  O  4J 

•r-  E 
+J  l- 

<D  E  O 

<U  +J  *r-  QJ 

-C  C  X  Q. 
+J  0)  co  co 
o  E  <u 

a>  s- 

C  <D  <D 
•«-  .C  -C  -C 
J.  P4J-P 

"O  *P4-  J 

(0  r— 

to  (0  c 

3  tO  -C  O 
(J  •!—  •!— 

O  4->  «4-> 
4-  co  ra  co 
S-  *r 
CO  3J=  S- 
E  U4J  (O 
to  O  “O  > 
CO  O  ••— 
r—  Jr- 
Q-r—  CO 
CO  f—  O 

a;  c  i—  *r- 

r  wop 

4-  to 
(O  *r 
C  «P  4-> 

i-  a;  co  co 

-XZ  4J 
CT>4->  C  tO 
C  O 
•i-  C  •«-  0) 
S-  t-  4J  -c 
S-  «3  4J 
3.X  S- 
U  CO  3  tO 
U  QJ  O  4J 
O  Q-  C 
<D  <D 
to  CD  C3  (O 
4J  -C  CO  QJ 
C  4J  S-  S- 
<U  CL)  Q- 

>  J=  >  QJ 
d)  O  fO  J— 

•r 

to  r  qj  a; 
3  ^JC  C 
O  4 J  •»- 
t-  4->  r— 
S-  CO  CO 
CO  4->  C 

>  0)  c  a; 
E  <u  -x 

4-  •!—  to  o 
o  +j  a;  u 

S-  X3 
to  C  Q. 

C  co  (D  (D 
o  o  i.r 
to  E  f- 

•I-  *o 

5-  Q)  •«- 

CO  x:  r-  • 
CXh—  O  r- 

E  to  CO 
O  C 

o  •  a;  a> 
x:  x:  •»- 

r-  UH  tO 
CO  C 

i-  -i-  a; 
o  a.  •  -c 
o.  a; 

E  fl)  c 

<D  -C  »r  4- 
h+Jr-  O 


R-band 


/ 


Figure  22.  Simultaneous  waveforms  comparing  the  time  of  the  R  band  emission 
with  respect  to  the  current  interruption  (!)  and  the  soft  x-ray 
emission.  The  initial  peaks  are  simultaneous  time  markers  whose 
positions  have  been  adjusted  to  reflect  differences  in  cable 


lengths  and  time-of-flight. 
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Figure  23. 
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Simultaneous  waveforms  comparing  the  occurrence  of  the  emission 
of  the  R  band  and  the  49  keV  portion  of  the  electron  beam  (60° 
channel  on  the  electron  spectrometer,  Figures  18  and  19).  The 
initial  peaks  on  the  waveforms  are  simultaneous  time  markers 
whose  positions  have  been  adjusted  to  reflect  differences  in 
cable  lengths  and  times-of-flight. 
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Figure  24.  Simultaneous  waveforms  comparing  the  occurrence  of  the  emission 
of  the  R  band,  the  83  keV  portion  of  the  electron  beam  (75° 
channel  on  the  electron  spectrometer.  Figures  18  and  19)  and  the 
hard  x  rays  (P.M.  tube).  The  initial  peaks  on  the  waveforms  are 
simultaneous  time  markers  whose  positions  have  adjusted  to  reflect 
differences  in  cable  lengths  and  times-of-f light.  The  last  peak 
on  the  PM  tube  signal  Is  due  to  2.45-MeV  neutrons  generated  by 
D(d,n)3He  reactions  in  the  device.  The  position  of  this  peak 
has  not  been  adjusted  for  time-of-flight. 


related  since  they  do  not  occur  at  the  sane  time.  This  correlation  might  have 

3 1 

been  expected  given  the  results  from  other  DPF  devices. 

B.2.4  Preliminary  Model 

« 

These  results  have  been  used  to  formulate  a  model  for  the  microwave 

emission  and  enhanced  resistivity  that  is  being  tested  at  present. 

The  essence  of  the  model  is  that  the  microwaves  and  enhanced  resistance  are 

caused  by  current  driven  instabilities  occuring  in  the  thin  current  sheath  (5 

plasma  skin  depths  thick)  produced  on  the  outer  boundary  of  the  pinch  during  1st 

compression.  This  is  similar  to  the  model  developed  by  Peacock  and  his  co- 

workers  to  explain  the  results  of  their  laser  scattering  and  Faraday  rotation 

experiments.  Their  results  indicate  that  at  the  maximum  in  the  1st 

compression  the  current  sheath  is  confined  to  a  thin  outer  skin.  The  current 

flow  redistributes  Itself  in  the  next  tens  of  nanoseconds  permitting  30%  of  the 

current  to  flow  in  the  interior.  Blue  shifts  of  the  scattered  light  indicate 

that  the  local  electron  drift  velocity  in  the  outer  skin  ranges  between  one  to 

one-tenth  times  their  thermal  speed.  Enhanced  scattering  from  this  region  indi- 
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cates  the  presence  of  ion  turbulence. 

This  model  for  the  experimental  results  from  Culham  is  further  supported  by 

laser  scattering  results  from  Frascati  on  their  last  1MJ  DPF  facility.  The 

larger  diameter  plasma  produced  by  that  device  (diameter  of  1.4  cm  versus  0.5  cm 

for  the  Culham  device)  permitted  the  Frascati  scientists  to  scatter  laser  light 

from  distinct  regions  to  the  plasma.  They  found  enhanced  scattering  from  the 

outer  regions  and  thermal  scattering  from  the  core  during  the  period  between 

3  8 

maximum  compressions  and  the  onset  of  neutron  production.  The  presence  of  a 
turbulent  outer  skin  current  following  first  compression  is  therefore  fairly 
well  established  experimentally. 


I 


Applying  the  model  to  our  results  the  following  predictions  will  be  compared 
with  experiment: 

1.  Frequency  of  the  microwave  emission 

2.  Resistance  of  the  pinch  versus  time 

3.  Duration  of  the  microwave  burst 

4.  The  timing  of  the  burst 

5.  The  power  level  of  the  burst 

The  results  of  these  comparisons  will  be  reported  next  fall  in  the  Annual 
Report  for  another  grant  (Grant  No.  AFOSR-84-0391 ). 
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ABSTRACT.  Particle  beam  diagnostic  tools  involving  an  ion  beam  Faraday  cup,  an  electron  beam  Faraday 
cup,  an  electron  magnetic  spectrometer  and  solid-stale  nuclear  track  detectors  have  been  used  to  measure  the 
parameters  of  the  particle  beams  generated  by  a  Mather-type  dense  plasma  focus.  At  a  capacitor  bank  energy 
of  12.5  kJ,  the  energy  spectra  of  the  electron  and  ion  beams  are  found  to  obey  the  same  power  laws: 
dN/dEocE"*  where  x  a  3.5  ±0.5.  Primary  electron  beam  current  and  energy  spectra  were  measured  as  a  func¬ 
tion  of  main  bank  current  ait  pinch  time,  Imb-  The  primary  electron  beam  current  was  found  to  scale  as 
•mb  aS>  reaching  a  magnitude  of  17  kA  for  a  device  energy  of  12.5  kJ.  The  exponent  x  of  the  electron  energy 
spectra  was  found  to  scale  as  x  <zrj^BSt0‘1.  These  results  ate  incorporated  into  an  axial  beam  target  model  for 
neutron  production,  and  it  is  found  that  this  model  could  account  for  the  magnitude  and  scaling  of  the  observed 
neutron  yield  with  Imb- 


1.  INTRODUCTION 

The  measurement  of  the  energy  spectra  of  particles 
accelerated  by  a  plasma  focus  (Fig.  1 )  is  of  interest  for 
several  reasons.  One  is  to  clarify  the  mechanism  of 
neutron  production  by  using  the  measured  ion  energy 
spectra  to  predict  the  beam-target  component  of  the 
neutron  yield.  A  second  reason  is  to  further  define  the 
acceleration  process  by  comparing  the  direction  and 
energy  spectra  of  the  various  particle  species  emitted. 
The  definition  of  the  acceleration  process  would  make 
possible  scaling  predictions  for  the  use  of  the  dense 
plasma  focus  (DPF)  as  an  accelerator  and/or  as  an 
inertial  confinement  fusion  driver.  Also  electrical 
energy  appears  to  be  coupled  to  the  plasma  focus  by 
current  interruption  ( I )  and  since  current  interruption 
has  been  proposed  as  mechanism  for  particle  acceler¬ 
ation  in  solar  flares  [2],  the  results  of  such  '•tudrs  may 
have  applications  to  astrophysics  as  well, 

*  This  article  Is  based  on  a  portion  of  the  PhD  thesis 
submitted  by  one  of  the  authors  (WS)  in  partial  fulfilment  of 
the  requirements  for  the  doctoral  degree  at  the  University  of 
Illinois. 

*  Present  address:  Sandia  National  Laboratory,  Albuquerque, 
New  Mexico,  87 185,  USA. 


IONS 


F1G.1.  Plasma  focus  device.  During  current  interruption  at 
pinch  time  the  ions  and  electrons  are  accelerated  in  directions 
shown. 


Figure  2  illustrates  that  current  interruption  is  a 
source  of  energy  in  a  plasma  focus.  The  lowest  trace 
is  the  main  bank  current  waveform  of  the  plasma 
focus.  The  sharp  reduction  in  the  current  by  a  factor 
of  two  at  pinch  time  indicates  [  1 1  that  one-half  the 
magnetic  energy  was  coupled  into  the  pinch  (i.e.  3.8  kJ). 
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FIG.2.  Simultaneous  waveforms  of  fast  electron  primary 
current  as  measured  on  fast  Faraday  cup  ( top ),  net  electron 
current  signal  as  measured  on  Rogowski  coil  mounted  on 
electron  drift  tube  (middle;  see  Fig.4)  and  current  flowing  in 
the  DPF  as  measured  by  a  Rogowski  coil  mounted  on  anode 
(bottom,  see  Fig.4).  Rectangular  pulses  at  100  ns  on  the  top 
two  traces  are  time  markers  and  electron  currents  are  simulta 
neous  with  end  of  dip  in  plasma  focus  current  (±50  ns).  Filling 
pressure  was  3  ton  deuterium;  the  bank  voltage  was  25  kV. 


FIGS.  Ion  Faraday  cup  time-of- flight  anangement.  Differen¬ 
tial  pumping  hole  is  0.48  cm  long. 


-OP F  HOLLOA-CENTRE 
ELECTRODE 


The  primary  current  of  the  electron  beam  generated 
during  the  DPF  current  interruption  is  also  shown  in 
Fig.  2.  This  is  the  signal  of  a  slightly  modified  Pellinen 
(3]  electron  beam  Faraday  cup  located  below  the 
device  (Fig.  1 ).  A  6-pm  mylar  filter  was  placed  over 
the  cup  to  permit  the  observation  of  the  primary  elec¬ 
tron  beam  current.  The  time  markers  show  that  the 
fast  electrons  are  generated  during  the  sharp  drop  in 
the  focus  current,  indicating  that  current  interruption 
is  converting  a  fraction  of  the  magnetic  energy  to 
particle  kinetic  energy. 

In  this  paper,  the  first  simultaneously  observed  energy 
spectra  of  ions  and  electrons  emitted  from  a  plasma 
focus  are  reported.  Also  the  scaling  of  the  primary 
electron  beam  current  and  energy  spectrum  with  the 
total  DPF  main  bank  current  at  pinch  time,  Imb>  *re 
reported  for  the  first  time  [4J.  The  ions  and  electrons 
are  accelerated  in  the  directions  shown  in  Fig.  1 .  The 
ion  beam  energy  spectra  were  measured  with  a  Faraday 
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FIG.  4.  Schematic  view  of  electron  magnetic  spectrometer. 
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cup  and  time-of-flight  analysis  [5]  (Fig.3)  and  the 
electron  beam  energy  spectra  with  a  magnetic  spectro¬ 
meter  (Fig.4). 

Previous  measurements  of  the  DPF  ion  beam  energy 
spectra  have  been  made  with  activation  analysis  [6—8], 
magnetic  spectrometers  (9,  10],  time-of-flight  analysis 
with  doubly  diffused  Si  (PIN)  detectors  [11,  12], 
Thomson  spectrometers  [13,  14],  and  nuclear  emulsions 
[15].  However,  the  time-of-flight  arrangement  used  in 
this  work  can  resolve  the  iorj  energy  spectrum  down 
to  lower  energies  [5]  than  any  of  the  other  methods, 
which  is  significant  since  the  electron  beam  energy 
spectrum  will  be  shown  to  have  a  large  component  in 
this  lower  energy  range. 

The  electron  beam  energy  spectrum  of  the  plasma 
focus  has  generally  been  inferred  by  analysis  of  the 
hard  X-ray  emission  spectrum  [  1 6—2 1  ]  or  of  dendrite 
patterns  [21-23]  caused  when  the  fast  electrons  strike 
a  target  such  as  plexiglas.  The  hard-X-ray  emission 
energy  spectra  are  observed  [16—19,21]  to  have  a 
power  law  dependence  decreasing  with  photon  energy, 
but  it  is  not  possible  to  uniquely  determine  the  electron 
energy  spectra  from  such  X-ray  measurements.  Usually, 
a  thick-target  bremsstrahlung  model  and  a  power-law 
electron  energy  spectrum  are  assumed  in  order  to 
evaluate  the  exponent  of  the  power  law  [16].  Analysis 
of  dendrite  patterns  in  target  materials  is  cumbersome 
and  the  lower  limit  of  the  observable  energy  is  about 
100  keV  [22]. 

The  magnetic  spectrometer  used  in  this  work  has 
none  of  these  drawbacks  and  while  this  type  of 
spectrometer  has  been  used  to  study  electron  beams 
generated  by  vacuum  diodes  [24,  25],  this  is  the  first 
time  it  has  been  used  to  study  the  fast  electrons  emitted 
by  a  plasma  focus. 

In  the  next  four  sections,  the  apparatus  and  the 
results  are  presented  and  discussed  and  conclusions 
relating  to  neutron  production  and  the  acceleration 
process  are  drawn. 

2.  THE  EXPERIMENTAL  APPARATUS 

The  Illinois  DPF  has  Mather-type  electrodes  as  shown 
in  Fig.l.  The  centre  electrode  (anode)  is  hollow  tb 
permit  direct  observation  of  fast  particles  emitted  from 
the  lower  end  of  the  pinch.  The  outer  diameter  of  the 
centre  electrode  is  5.08  cm  and  the  inner  diameter  of 
the  squirrel  cage  (outer  electrode)  is  10.2  cm.  The 
centre  electrode  extends  14  cm  beyond  a  lavite 
insulator  [26].  The  electrode  material  is  OFHC  copper; 
there  is  evidence  that  oxygen  released  from  regular 


copper  during  a  discharge  is  detrimental  to  focus 
formation  [27],  The  DPF  capacitor  bank  (40  pF, 

60  kV)  is  composed  of  four  modules,  each  of  which  is 
switched  by  a  trigatron  spark  gap.  A  trigger  pulse 
generator  has  been  built  [28]  which  provides  80  kV 
pulses  with  a  1 0-90%  rise-time  of  less  than  5  ns  to 
trigger  the  switches.  The  jitter  in  the  switch  closure 
time  is  measured  to  be  less  than  4  ns.  The  capacitor 
bank  modules  are  coupled  to  the  DPF  electrodes  via 
four  low-inductance  parallel-plate  transmission  lines 
arranged  symmetrically  around  the  focus.  The  total 
circuit  inductance  external  to  the  DPF  electrodes  is 
1 7  nH.  The  DPF  chamber  is  pumped  down  to  1  fin l 
between  shots  by  a  trapped  two-stage  roughing  pump. 

For  the  measurements  reported  here  the  bank 
voltage  was  16  to  25  kV,  the  deuterium  fill  gas  pressure 
was  3  torr  and  the  bank  current  at  pinch  time  was  up 
to  560  kA.  Under  these  conditions  the  average 
neutron  yield  is  3.8  X 109  and  the  peak  yield  is 
7.5  X 10*.  At  constant  pressure  (3  torr  deuterium), 
the  neutron  yield  is  found  to  scale  as  the  plasma  focus 
current  at  pinch  time  to  the  power  of  4.4  (Fig.5).  The 
bank  voltage  was  limited  to  25  kV  and  below  for  the 
results  reported  here.  This  limit  was  selected  because 
fiashover  damage  to  the  capacitor  bank  insulation 
occurred  at  this  voltage  when  the  device  was  operated 
at  a  fill  pressure  of  1  torr  D2 .  The  device  was  normally 
operated  at  a  higher  pressure  (3  torr)  which  is  predicted 


(LOG  SCALE) 


FIG.5.  Plot  of  neutron  yield  versus  plasma  focus  current  at 
pinch  time.  Filling  pressure  was  5  torr  deuterium;  bank  voltage 
range  16-24  kV. 
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[29]  to  lower  the  induced  voltage  across  the  device  and 
hence  provide  some  margin  of  safety. 

The  energy  spectra  of  the  ions  and  electrons  reported 
here  are  those  of  particles  directed  along  the  symmetry 
axis  of  the  device  and  in  the  direction  of  the  current 
flow  for  each  species.  Since  the  outer  electrode  is 
negative  with  respect  to  the  centre  electrode  the  ions 
move  away  from  the  centre  electrode  and  the  electrons 
pass  through  it  (Fig.  1 ).  The  ion  energy  spectrum  was 
measured  with  a  magnetized  Faraday  cup  and  time- 
of-flight  analysis;  details  are  reported  elsewhere  [5  J. 
The  basic  feature  of  the  time-of-flight  arrangement  is 
the  use  of  a  re-entrant  drift  tube  with  a  0.046-cm- 
diameter  hole  on  the  end  placed  1 6  cm  above  the  anode 
on  the  device  axis  (Fig. 3).  Most  of  the  48  cm  ion  flight 
path  is  differentially  pumped  to  1 5  jam  and  thus  the 
effects  of  energy  attenuation  and  multiple  scattering 
of  the  ions  by  the  gas  molecules  are  greatly  reduced, 
allowing  resolution  of  the  ion  energy  spectra  down  to 
25  keV. 

The  instant  of  acceleration  for  time-of-flight  analysis 
is  taken  to  be  the  onset  of  the  hard  X-rays  observed 
with  a  plastic  scintillator-photomultiplier  tube  com¬ 
bination.  It  is  assumed  that  the  electrons  accelerate 
at  the  same  instant  as  the  ions  and  produce  thick- target 
bremsstrahlung  when  striking  the  centre  electrode  or 
the  vacuum  chamber.  Although  the  analysis  assumes 
that  only  one  ion  species  (i.e.  D*)  is  present,  this 
assumption  is  consistent  with  the  observation  that 
with  deuterium  gas  only  one  peak  in  the  Faraday  cup 
signal  is  observed  [5]  and  nearly  equal  amounts  of 
hydrogen  and  deuterium  must  be  present  for  two 
peaks  to  be  observed  [5],  In  the  latter  case,  the  ratio 
of  the  flight  times  of  the  two  peaks  is  ~  1.4  which  is 
consistent  with  the  species  H*  and  D*  undergoing  ion 
diode  acceleration  [30]  followed  by  multiple  scattering 
in  the  gas  target  ahead  of  the  beam  (see  Section  4). 

To  interpret  the  ion  Faraday  cup  signals,  corrections 
are  made  [5]  for  charge  state,  multiple  scattering  by 
the  gas  and  energy  attenuation. 

To  measure  the  electron  beam  energy  spectrum,  the 
electron  beam  is  transported  from  the  DPF  to  the 
magnetic  spectrometer  as  shown  in  Fig. 4.  The  drift 
space  diameter  of  the  25  cm  long  hollow  anode  is 
3.18  cm;  the  diameter  of  the  45  cm  long  aluminium 
drift  tube  is  4.45  cm.  Similar  dimensions  have  been 
used  previously  to  transport  electron  beams  (with 
vly  <  1  and  v/y  =»  1)  with  less  than  20%  loss  of  beam 
charge  over  distances  of  ~  1  m  at  optimum  ,  ressure 
[31,  32].  The  optimum  pressure  is  defined  to  be  that 
pressure  at  which  the  electron  beam  charge  is  trans¬ 
ported  with  minimum  loss  [33,  34],  and  is  also  the 


pressure  at  which  a  beam  is  transported  with  the  least 
amount  of  change  in  its  energy  spectrum  [35].  For 
deuterium,  the  optimum  pressure  range  (80%  trans¬ 
mission)  for  beams  with  vh  *  1  is  from  2  to  5  torr 
[33,  34];  consequently,  the  same  gas  used  in  the 
DPF  chamber  (i.e.  3  torr  deuterium)  is  also  used  in  the 
drift  regions. 

The  characteristics  of  the  magnetic  spectrometer 
are  as  follows: 

1 .  The  drift  tube  and  spectrometer  chamber  are 
made  from  1/4  in.  thick  aluminium  to  shield 
against  external  fields  and  to  minimize 
scattering  and  X-ray  production. 

2.  Tantalum  collimators  limit  the  collected  current 
to  a  few  milliamperes  to  prevent  collective 
effects  from  disrupting  electron  orbits  inside  the 
spectrometer. 

3.  An  Armco  shield  reduces  the  magnet’s  fringe 
field  in  the  collimating  region  by  a  factor  of  10-3 

4.  1 0  spectrometer  channels  are  available  that 
cover  energies  from  ~  10  keV  to  ~10  MeV  for 
a  200-G  field. 

5.  Faraday  cup  detectors  are  used;  consequently, 
time  resolution  is  limited  by  the  bandwidth  of 
the  oscilloscope  channels  (~100  MHz). 

6.  Time-integrated  spectra  can  also  be  obtained  by 
appropriate  use  of  the  signal  cable  capacitance. 

7.  The  Faraday  cup  collectors  are  situated  in  the 
fringe  field  of  the  electromagnet  to  inhibit 
secondary-electron  emission. 

The  spectrometer  was  calibrated  by  mapping  the 
magnetic  field  produced  by  the  electromagnet  with  a 
gaussmeter  and  using  a  computer  code  to  follow 
particle  trajectories  inside  the  spectrometer.  The  error 
due  to  the  finite  width  of  the  collimators  is  calculated 
to  be  about  1 0%. 

An  electron  beam  Faraday  cup  was  used  to  obtain 
a  direct  quantitative  measurement  of  the  primary 
current  of  the  electron  beam  generated  by  the  plasma 
focus.  Rogowski  coils  have  been  used  to  estimate  the 
magnitude  of  the  electroi  beam  current  [36-39],  but, 
because  of  the  plasma  return  current,  interpretation 
of  Rogowski  coil  signals  is  not  straightforward. 

Figure  2  shows  simultaneous  wave  forms  of  Rogowski 
coil  and  Faraday  cup  waveforms.  The  Faraday  cup  has 
a  thin  mylar  window  to  inhibit  the  flow  of  plasma 
return  current. 

The  electron  beam  Faraday  cup  that  was  built  is  a 
slightly  modified  [4]  Pe.linen  [3]  Faraday  cup.  The 
cup  was  designed  to  measure  1 00-kA  currents,  has  an 
experimentally  determined  [4]  bandwidth  of  450  MHz, 
and  the  1  /4  mil  mylar  plasma  current  filter  normally 
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used  permits  primary  electrons  with  energies  down  to 
27  keV  to  be  collected  by  the  cup. 

A  solid-state  nuclear  track  detector  (SSNTD)  tech¬ 
nique  to  measure  high  fluences  of  light  ions  [40]  was 
developed  as  a  means  to  quantitatively  measure  the 
total  number  of  ions  accelerated.  Fast  deuterons 
striking  a  converter  layer  of  Li2B407  interact  to  pro¬ 
duce  energetic  (d,  a)  alpha  particles  which  are  recorded 
by  a  SSNTD  placed  on  the  opposite  side  of  the  layer. 

If  the  shape  of  the  deuteron  energy  spectrum  is  known, 
the  track  density  can  be  used  to  estimate  the  incident 
beam  fluence.  The  ion  energy  spectral  shape  deter¬ 
mined  by  the  ion  Faraday  cup  and  the  observed  track 
density  indicate  that  1  kJ  of  the  energy  deposited  in 
the  plasma  (~3.8  kJ  out  of  12.5  kJ)  goes  into  deuterons 
accelerated  to  25  keV  or  more.  The  corresponding 
energy  efficiency  for  this  process  is  about  8%. 

The  four  diagnostics  described  enable  both  a  qualita¬ 
tive  estimate  of  the  electron  and  ion  energy  spectra 
and  a  quantitative  estimate  of  the  total  fluence  of  the 
beam  particles.  The  results  of  the  experiments  using 
these  diagnostics  are  reported  in  the  next  section. 


E«(k«V) 


FIG.  7.  Time-integrated  electron  beam  energy  spectra  for  four 
separate  plasma  focus  shots  taken  on  3  Nov.  1980.  Data  were 
normalized  to  lowest  energy  channel  (47  to  74  keV).  Charging 
voltage  was  25  kV  (12.5  kJ);  pressure  of  deuterium  fill  gas  was 
3. 0  ton,  and  neutron  yield  for  all  these  shots  was  within  20% 
of  3.8X10*. 
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FIG.  6.  Ion  energy  spectra  from  four  shots  taken  on  25  June 
1980.  Charging  voltage  on  the  capacitor  bank  was  24  kV 
(11.5  kJ);  pressure  of  deuterium  fill  gas  was  3.0  ton,  and 
neutron  yield  for  all  these  shots  was  within  10%  of  2.2  X 10*. 
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FIG. 8.  Simultaneously  recorded  ion  and  electron  energy 
spectra.  Straight  lines  are  power  law  fits  to  the  ion  and  electron 
spectra  where  exponents  are -3.5  and -3.8,  respectively. 


3.  RESULTS 

The  results  of  the  energy  spectra  measurements  are 
shown  in  Figs  6-8.  In  Fig.6,  the  ion  spectral  data  for 
four  DPF  shots  are  shown  for  a  deuterium  fill  gas 
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pressure  of  3  torr  and  a  bank  voltage  of  24  kV.  The 
neutron  yields  for  these  four  shots  are  within  1 0%  of 
2.2  X 1 0* .  Although  there  is  some  evidence  of  a  bump 
on  the  tail  the  overall  power  law  spectra  are  about  the 
same.  The  exponent  of  this  power  law  is  -3.4  ±  0.4. 

The  electron  energy  spectra  were  measured  with  the 
magnetic  spectrometer  for  four  shots  under  similar 
conditions  (3  torr  deuterium,  25  kV  bank  voltage). 

The  spectrometer  was  operated  in  the  time-integrated 
mode  and  the  results  are  shown  in  Fig.7.  As  with  the 
ions,  power-law  energy  spectra  are  observed  on  each 
of  the  shots  with  a  mean  value  for  the  exponent. of 
-3.6  ±  0.3.  The  neutron  yields  for  these  four  shots 
lie  within  20%  of  3.8  X  10*.  This  power  law  is  equi¬ 
valent  to  that  inferred  by  Lee  et  al.  [  1 6]  and 
Van  Paasen  et  al.  [  18]  from  the  hard  X-ray  spectra 
emitted  from  plasma  foci  of  similar  energies  and 
neutron  yields.  It  is  also  interesting  that  this  same 
spectrum  of  hard  X-rays  was  observed  during  the 
reconnection  phase  of  a  double  inverse  pinch  [41  ], 
possibly  indicating  a  similar  process  in  a  device  of  quite 
different  geometry. 

The  results  of  a  simultaneous  observation  of  the 
electron  and  ion  beam  energy  spectra  are  shown  in 
Fig.8.  The  conditions  are  the  same  as  described  in 
the  preceding  paragraph  and  the  neutron  yield  was 
3.4  X  1 09 .  The  two  spectra  have  the  same  power  law 
structure  to  within  the  uncertainty  in  the  measure¬ 
ments  and  are  in  agreement  with  the  results  of  the 
other  measurements  at  these  conditions  shown  in 
Figs  6  and  7. 

The  electron  energy  spectrum  and  total  primary 
electron  current  were  observed  for  main  bank 
currents  between  350  and  560  kA.  To  avoid  the 
complications  of  radial  spoke  formation  and  parasitic 
currents  [42]  which  may  vary  from  shot  to  shot  only 
those  data  were  selected  which  had  a  neutron  yield 
within  a  factor  of  two  of  the  best  yield  at  that  bank 
voltage.  This  occurs  over  60%  of  the  time  at  25  kV 
where  the  spark  gaps  and  anode  length  are  more 
optimal  but  occurs  less  frequently  when  the  bank 
voltage  is  reduced.  It  is  assumed  that  the  effects  of 
parasitic  current  formation  are  minimal  at  high  neutron 
yields  (there  is  no  neutron  yield  when  radial  spokes 
form  [42])  and  the  peak  yield  scales  as  the  I^1 0,3  for 
the  Illinois  plasma  focus  (Fig.5),  a  dependence  reported 
by  several  other  investigators  when  the  performance  of 
their  plasma  focus  devices  is  optimal  [43]. 

The  electron  magnetic  spectrometer  was  operated  as 
shown  in  Fig.4.  The  perpendicular  magnetic  field  was 
nominally  at  200  G  and  six  energy  channels  were  used, 
spanning  a  mean  energy  range  between  30  and  400  keV, 


FIG.  9.  Typical  time-integrated  electron  energy  spectrum 
observed  with  electron  magnetic  spectrometer  (see  Fig.4)  for 
IM-S60kA.  Filling  pressure  was  3  torr  deuterium;  bank 
voltage  was  25  kV  ((1/2)  CV1  - 12.5  kl). 


FIG.  1 0.  Typical  time-integrated  electron  energy  spectrum 
observed  with  the  electron  magnetic  spectrometer  (tee  Fig.4) 
forlufo*370kA.  Filling  pressure  was  3  torr  deuterium:  bank 
voltage  was  17.5  kV  ((1/2)  CK1  -  6. 1  kj). 
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OPF  MAN  SANK  CURRENT  AT  PINCH  TIME 
IN  KILOAMPS  (LOG  SCALE) 

FIG.  11.  Observed  scaling  of  the  power  law  exponent, 
x(dN/dE<xE'%  ),  with  /mb-  Filling  pressure  was  3  torr  deuterium; 
vertical  error  bars  represent  least-square  error  in  power  law 
slopes  for  various  observations 


As  above,  power  law  energy  spectra  were  observed 
for  the  electrons  (Fig.9)  although  bumps-on-the-tail  of 
the  energy  distributions  were  sometimes  observed  at 
lower  bank  currents  (Fig.  1 0).  Since  the  bump  forma¬ 
tion  process  appears  to  be  different  from  the  power 
law  formation  process  (44]  and  is  intermittent,  these 
bumps  were  not  included  in  calculating  power  law 
exponents.  The  results  are  shown  in  Fig.  1 1 .  The  elec¬ 
tron  spectra  appear  to  harden  with  increasing  I  mb 
where  the  absolute  value  of  the  exponent,  |x|,  scales 
as  Imb’±#  *.  These  results  are  consistent  with  the 
hard-X-ray  data  where  the  absolute  value  of  the  electron 
beam  energy  power  law  ranged  from  3.5-4.0  for 
20  kJ  devices  (16,  18]  to  2.4-4.0  for  a  57.3  kJ  device 
(17],  and  down  to  2.0  for  a  375  kJ  device  (19].  Thus 
the  lesults  from  the  6  - 12.5  kJ  device  reported  here 
fit  these  X-ray  data  fairly  well  although  these  new 
results  are  less  ambiguous,  being  the  result  of  a  direct 
measurement. 

The  fill  pressure  for  the  DPF  for  all  the  data  reported 
here  is  3  torr  D, .  However,  varying  the  pressure 
between  2  and  4  torr  at  a  bank  voltage  of  25  kV 
(1mb~  560  kA)  had  no  effect  on  the  power  law 
(-3.4  i  0.4)  although  the  signals  were  stronger  at 
2  torr  and  weaker  at  4  torr. 

The  scaling  of  the  peak  in  the  primary  current  with 
Imb  was  also  determined.  The  electron  magnetic 
spectrometer  (Fig.4)  was  replaced  by  the  electron 
Faraday  cup  (Fig.  12)  to  measure  the  primary  current 
of  the  electron  beam.  The  filter  had  to  be  replaced 
after  each  shot  and  so  a  gate  valve  was  inserted  into 
the  drift  tube  and  opened  after  the  plasma  focus  had 
been  conditioned. 
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FIG.  12.  Schematic  view  of  high-current  fast  Faraday  cup  used 
to  measure  primary  electron  current.  Areas  behind  filter  are 
evacuated  to  about  0.04  torr  to  increase  voltage  stand-off. 


DPF  MAIN  BANK  CURRENT  AT  PINCH  TIME 
IN  KILOAMPS  (LOG  SCALE) 

FIG.  1 3.  Observed  scaling  of  peak  primary  current  (first  peak 
In  top  trace  of  Fig. 2)  with  /mb-  Filling  pressure  was  3  torr 
deuterium. 


Typical  data  are  shown  in  Fig.2  where  the  Rogowski 
coil  signal  (net  current)  is  included  for  reference.  When 
the  filter  is  removed,  the  Faraday  cup  signal  closely 
follows  the  form  of  the  Rogowski  coil  signal  although 
reduced  by  about  15%,  presumably  because  of  beam 
divergence.  The  difference  between  these  two  signals 
illustrates  the  difficulty  in  interpreting  Rogowski  coil 
signals  caused  by  these  electron  beams  which  have 
been  reported  in  the  past  (36-39]. 

The  results  of  the  scaling  measurements  for  the 
primary  beam  current  are  shown  in  Fig.  13.  The  data 
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indicate  that  the  primary  beam  current  increases  with 
I  MB  to  a  power  of  2.9  ±0.5.  These  measurements  are 
the  first  of  their  kind  although  the  scaling  of  the  net 
current  with  Imb  has  been  previously  reported  by  this 
group  with  a  similar  power  law  (39).  The  data  point 
reported  by  Molen  (2 1  ]  is  included  for  reference. 

An  estimate  of  the  total  number  of  electrons  in  the 
fast  peak  observed  on  the  Faraday  cup  at  12.5  kJ  is  of 
order  2  X 1 0l! .  The  mean  ion  beam  fluence  over  a 
1 5-cm-diameter  plate  20  cm  away  from  the  focus 
estimated  using  the  SSNTD  technique  [40]  for  a  similar 
shot  is  6  ±  2  X  10,4cm"* .  Since  the  track  density  is 
observed  to  be  uniform,  we  obtain 

j-7-  ~  50  ±  30 

Ne 

where  N,  and  Ne  are  the  numbers  of  ions  and  electrons 
accelerated,  respectively.  The  ratio  is  of  order 
(M^Me)17* ,  which  js  predicted  by  a  model  for  a 
uniform-current-density  diode  accelerator  [45]. 

The  assumptions  that  fast  deuterons  are  created  at 
the  same  instant  as  the  fast  electrons  (or  the  X-ray 
bremsstrahlung)  and  at  the  top  of  the  anode  are  consis¬ 
tent  with  three  experimental  observations  [5]. 

1.  The  ratios  of  the  flight  times  for  the  two  peaks 
observed  when  mixtures  of  deuterium  and 
hydrogen  are  used  as  a  factor  of  yfl. 

2.  The  peaks  in  the  raw  time-of-flight  data  occur 
at  80  keV  as  consistent  with  multiple-scattering 
theory  [46]  in  3  torr  D2  with  a  path  length  of 
6  cm  which  is  inferred  from  geometrical 
considerations. 

3.  The  measured  secondary-emission  coefficient 
for  copper  has  the  same  energy  dependence  as 
that  observed  in  dedicated  experiments  [47]. 

Since  these  observations  involve  independent  pheno¬ 
mena  the  self-consistency  is  fairly  solid  evidence  for 
the  simultaneous  creation  of  the  two  species  within 
experimental  uncertainties  (±2cm,  ±5  ns). 

Figures  14  and  15  show  that  the  fast  electrons  as 
evidenced  by  the  100-keV  channel  of  the  electron 
magnetic  spectrometer  and  the  hard-X-ray  brems¬ 
strahlung  peak  on  our  plastic  scintillator-photomulti¬ 
plier  tube  occur  simultaneously  with  a  peak  in  |1| .  The 
bremsstrahlung  peak  occurs  at  the  second  and  slightly 
smaller  |I|  peak  (Fig. 1 5)  and  this  may  indicate  a  change 
in  sign  in  L[29]. 

Figure  1 6  shows  the  magnitude  of  the  terms  not 
involving  L(V  =  q/C  +  Li)  estimated  for  Imb  =  560  kA. 
To  compare  this  with  the  maximum  (Em  )  energy  of  the 
particles  accelerated,  the  maximum-energy  channel 
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FIG.  14.  Simultaneous  waveforms  of  a)  100-keV  (±20  keV) 
channel  of  electron  magnetic  spectrometer  and  b)  1-waveform 
of  the  mainbank  current.  Filling  pressure  was  3  torr  deuterium; 
bank  voltage  was  25  kV. 
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FIG.  15.  Simultaneous  waveforms  of  a )  scin  dilator  photo¬ 
multiplier  tube  and  b)  /  waveform  of  main  bank  current.  Long 
narrow  peak  in  a)  has  been  identified  as  hard-X-ray  brems¬ 
strahlung  produced  by  the  fast  electrons  in  the  beam  and  the 
broader  second  peak  as  the  neutrons  produced  by  the  device. 
Waveform  a )  is  delayed  40  ±4  ns  with  respect  to  b),  due  to 
Inherent  delay  of  photomultiplier  tube,  differences  in  cable 
length  and  time-of-fiight.  Filling  pressure  was  3  torr  deuterium; 
bank  voltage  was  25  kV. 
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FIG.  1 6.  Typical  waveforms  of  a)  voltage  across  the  capacitor 
bank  (q/C)  and  b)  I  of  the  mainbank  current  (the  two  peaks 
are  simultaneous).  These  waveforms  can  be  used  to  calculate 
two  of  the  three  terms  of  the  possible  accelerating  voltage  across 
the  plasma  diode  (excluding  the  LI  term).  The  circuit  induc¬ 
tance  at  pinch  time  Is  estimated  to  be  50  nH  on  the  basis  of 
time-integrated  pin-hole  camera  pictures  so  the  total  of  q/C  and 
U  is  IIO  kV.  (The  base  line  in  a)  is  one  division  below  the  top. ) 


receiving  signal  is  the  400-keV  channel.  Thus,  the 
1 10  keV  estimated  in  Fig.  16  is  only  about  one  fourth 
of  Em  .  Evidence  from  measurements  using  solid-state 
nuclear  track  detectors  [40]  indicates  that  some 
deuterons  acquire  an  energy  greater  fb«n  1 .0  MeV 
Particles  with  energies  from  1-5  Me\  Lave  been 
observed  in  several  other  DPF  devices  [6—8,  11  —  15]. 


4.  DISCUSSION 

The  results  reported  in  Section  3  have  implications 
on  two  aspects  of  DPF  phenomena: 

1.  the  neutron  production  mechanism,  and 

2.  the  acceleration  mechanism. 

These  aspects  will  be  treated  in  this  section. 

Since  the  total  number  of  particles  accelerated  and 
their  energy  spectra  have  been  determined  one  can 
use  this  information  to  predict  the  magnitude  and  the 
scaling  of  the  axial  beam  target  neutron  yield  with  Imb 
if  a  few  simple  assumptions  are  made  about  the  target. 
If  the  target  is  assumed  to  be  the  cold  gas  ahead  of  the 
beam  then  it  is  independent  of  bank  parameters. 


TABLE  I.  THEORETICAL  AND  EXPERIMENTAL 
DATA 


Phenomena 

Experiment 

Cold  large! 

Hoi  largel 

V,«1MB.X 

4.6  ±0.3 

5.2  ±1.0 

4.8  ±1.0 

Targel  Ihickness 

15  cm,  3  lorr 

2.9  X  1 0,adeuterons*  cm”a 

Y„  at  12.5  kj 

3.8X10* 

2.4  X  10* 

3.1X10* 

Anisolropy 

1.5 

2.2 

1.7 

<Y,o  -/V.*-) 

Neutron  pulse 

90  ns 

60  ns* 

30  ns* 

width  at  3.2  m 

Assumes  axial  beam  motion  only  and  includes  20-ns  pulse  width 


DPF  MAW  BAMS  CURRENT  AT  PINCH  TIME  IN  KILOARPS 
(LOG  SCALE) 

FIG.  1 7.  Neutron  yield  In  Illinois  plasma  focus  and  predictions 
of  various  models: 

I)  Curve  a)  is  least-square  fit  to  peak  neutron  yield  versus 
^MB  ^  ^n  a  f  MB  *  °  *  A 

il)  Curve  b)  is  least-square  fit  to  neutron  yield  for  shots 
actually  used  to  determine  the  data  in  Figs  II  and  13. 
Triangles  are  data  with  representative  error  bars  shown  for 
one  of  the  points.  For  these  data  Yn®/jj,'010-3. 
iff  Curve  c)  is  result  of  beam  hot  target  model  for  a  target 
thickness  of  2.9  X  / 01*  deuterons- cm'1 .  Error  bars  are 
those  associated  with  experimental  uncertainties  in  the 
beam  parameters  shown  in  Figs  1 1  and  13.  The  model . 
predicts  Yn<x / **. 

iv)  Curve  d)  >s  result  of  beam  cold  gas  target  for  a  target  of 
3  ton  deuterium  gas  15  cm  long  (target  thickness 
2.9  X  10 '*  deuterons- cm'1  as  in  iii)  where  effects  of 
slowing-down  on  deuteron  beam  have  been  Included. 

Enor  bars  are  as  In  HI)  and  this  latter  model  predicts 


NUCLEAR  FUSION,  Vol.22.  No.9  ( 19*21 


1.169 


vv  •«: v 


Similarly,  if  one  assumes  a  hot  target  (plasma)  of 
density  and  length  independent  of  Imb.  which  is 
roughly  consistent  with  several  experiments  [48—52], 
even  slowing-down  of  the  deuterons  could  be  neglected 
for  Te>500  eV  for  targets  less  than  10  cm  in  length. 
Thus  in  calculating  the  axial-beam  target  neutron  yield 
the  following  assumptions  have  been  made: 

1 .  The  ion  beam  energy  spectrum  has  the  same 
power  law  dependence  as  the  electron  beam 
energy  spectrum  for  a  given  ImB' 

2.  N+/Ne  =  (MJMe)1/1. 

3.  The  ion  and  electron  energy  beam  spectra  are 
the  same  at  all  angles  to  the  axis  where  the  beams 
exist. 

4.  The  target  parameters  are  independent  of  Imb  • 

5.  The  slowing-down  in  the  cold  gas  target  is  that 
for  deuterons  in  3  torr  D3  gas  [53J. 

6.  The  slowing-down  in  the  hot  plasma  target  is 
negligible. 

The  predictions  of  this  model  are  compared  with 
the  experimental  data  in  Table  I. 

The  scaling  and  magnitude  of  the  beam  target 
neutron  yield  is  consistent  with  either  target  model 
(Fig.  17)  and  the  electron  beam  scaling  parameters 
reported  in  Section  3.  However,  neither  target  model 
is  adequate  to  explain  either  the  observed  anisotropy 
or  the  neutron  pulse  width.  The  long  pulse  width  and 
low  anisotropy  observed  experimentally  could  be  due 
to  trapping  of  part  of  the  beam  in  some  magnetic 
structure  of  the  focus  as  evidenced  by  the  high  ratio  of 
DT  neutrons  to  DD  neutrons  observed  in  other  DPF 
experiments  [50,  54],  Description  of  such  trapping  is 
beyond  the  scope  of  these  target  models. 

Owing  to  the  assumptions  made  in  the  calculations, 
the  agreement  in  the  magnitude  and  scaling  of  the  beam 
target  model  prediction  of  neutron  yield  and  the 
experimental  result  only  means  that  the  axial  beam 
target  model  cannot  be  ruled  out.  Probably  the  most 
questionable  assumption  is  that  the  deuteron  beam 
energy  spectrum  is  the  same  at  all  angles  (to  1 5°)  with 
respect  to  the  device  axis.  This  is  apparently  true  for 
the  electron  beam  which  is  not  always  on  axis  (pre¬ 
sumably  because  of  the  frozen-hose  instability  [55]) 
and  yet  always  gives  the  same  power  law  spectrum 
(within  the  error).  The  large  uniform  patterns  of 
tracks  observed  in  the  SSNTD  converter  layer  tech¬ 
nique  [40]  (which  has  a  lower  energy  threshold  of 
about  150  keV)  over  a  19°  angle  with  respect  to  the 
axis  also  support  the  assumption.  And,  finally,  the 
deuteron  energy  spectra  are  roughly  independent  of 
the  angle  for  a  28-kJ  DPF  at  5  torr  Da  in  the  energy 
range  from  0.9  to  4.5  MeV  [52]. 


One  possible  paradox  is  that  the  instantaneous  ion 
current  is  estimated  to  be  eN^/At  ~  800  kA  (where  the 
pulse  width  is  taken  to  be  20  ns)  and  the  peak  device 
current  is  only  5.6  X  10s  A.  However,  the  mean  ion 
speed  (2.1  X  10*cm*s-1 )  corresponds  to  an  electron 
energy  of  ~  1 3  e  V,  so  presumably  the  ions  drag  along 
low-energy  electrons  which  cancel  any  net  current  and 
corresponding  magnetic  energy. 


ACCELERATION  MECHANISM 

It  is  clear  from  Figs  14  and  15  that  the  electron  (and 
hence  ion)  beam  occui  at  the  peak  in  II  mb  I  >  indicating 
that  current  interruption  plays  a  role  in  the  acceleration 
process.  However,  estimates  of  the  induced  voltage 
across  the  pinch  from  the  observed  voltage  surge  across 
the  parallel-plate  transmission  lines  between  the  capa¬ 
citor  banks  and  the  DPF  device  and  LI  (where  L  is 
estimated  from  pinhole  camera  pictures)  indicate  that 
only  about  100  kV  appears  across  the  pinch  due  to 
IImbI  phenomena.  This  voltage  is  much  too  low  to 
explain  the  high  particle  energies  observed  by  simple 
diode  acceleration  if  current  interruption  (I)  was  the 
only  source  of  the  applied  potential. 

It  has  been  suggested  [56]  that  motion  of  the  current 
sheath  could  cause  induced  voltages  much  higher  than 
LI  effects  and  the  computational  results  of  Kondoh  and 
Hirano  [29]  contain  estimates  of  the  magnitude  of  this 
effect.  The  system  of  equations  used  by  the  latter 
authors  has  been  programmed  and  used  to  estimate 
the  induced  voltage  for  the  Illinois  DPF  device.  The 
final  radial  collapse  (r0  =  2.5  cm)  was  modelled  and  the 
initial  conditions  were  determined  experimentally  by 
placing  a  co-axial  short  across  the  open  end  of  the 
plasma  focus  electrode  and  measuring  the  circuit 
characteristics  when  the  device  was  evacuated.  At 
3  torr  Dj  pressure,  the  code  estimates  total  induced 
voltages  range  up  to  ~  100  kV  for  a  peak  device 
current  of  560  kA.  Hence  it  appears  that  even  sheath 
motion  cannot  provide  sufficient  potentials  to  explain 
the  observed  particle  energies  (which  apparently 
range  above  1  MeV  in  this  device  [40])  by  the  simple 
process  of  diode  acceleration. 


CONCLUSIONS 

The  energy  spectra  of  the  fast  ions  and  electrons 
emitted  by  the  plasma  focus  in  opposite  directions 
(Fig.  1 )  are  observed  by  direct  methods  to  have  the 
same  power  law  energy  dependence  at  a  device  current 
at  pinch  time,  Imb.  of  560  kA.  The  scaling  of  the 
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absolute  value  of  the  exponent  of  the  power  law  was 
found  (for  electrons)  to  decrease  with  increasing  I  mb 
(hardening  of  the  spectrum)  and  the  values  were  found 
to  be  consistent  with  the  electron  energy  spectra  inferred 
from  hard  X-ray  spectral  measurements  on  several 
other  devices  (16-19],  It  is  interesting  that  power 
law  spectra  of  similar  exponents  have  also  been 
observed  in  a  double  inverse  pinch  (41],  cosmic  rays 
(57],  and  solar  flares  (58], 

The  magnitude  of  the  fast  .Mectron  current  has  been 
measured  by  a  filtered  fast  Fa.  a  Jay  cup  and  found  to 
scale  as  Imb  to  the  third  power.  The  highest  primary 
beam  current  observed  was  1 7  kA  for  a  bank  energy  of 
12.5  kJ.  Since  much  more  energetic  plasma  foci  exist 
such  as  the  1  MJ  device  at  Frascati  (50],  it  would  be 
very  interesting  if  similar  measurements  were  performed 
on  these  devices  to  see  if  the  scaling  observed  in  the 
Illinois  device  persists  at  high  bank  energy.  If  so,  the 
hardening  of  the  energy  spectra  indicated  by  these 
measurements  and  others  (16-19]  may  lead  to  new 
applications  of  the  DPF  as  a  pulsed  electron  beam 
source. 

The  measured  beam  parameters  and  scaling  laws  can 
be  used  to  predict  the  scaling  and  magnitude  of  the 
neutron  yield  due  to  a  beam  target  model.  If  the 
target  parameters  are  assumed  to  be  independent  of 
Imb,  then  the  agreement  between  beam  target  yield 
and  the  observed  yield  is  quite  good  with  respect  to 
scaling  and  magnitude  but  the  predictions  of  the  model 
do  not  fit  the  observed  temporal  neutron  pulse  width 
and  anisotropy.  It  is  possible  that  trapping  of  all  or 
part  of  the  deuteron  beam  by  magnetic  structures  in 
the  plasma  focus  (50]  could  explain  the  latter 
discrepancies. 

The  beams  are  observed  to  occur  at  the  peak  in  the 
current  interruption  (Imb),  indicating  that  it  has  a 
strong  influence  on  the  acceleration  process  but  the 
magnitudes  of  the  voltages  generated  by  d/dt  (LI) 
effects  estimated  from  circuit  measurements  are  too 
low  to  explain  the  acceleration  process  by  simple  diode 
acceleration. 
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have  applications  to  astrophysics  as  well. 


f  This  article  is  based  on  a  portion  of  the  PhD  thesis 
submitted  by  one  of  the  authors  (WS)  in  pirtial  fulfilment  of 
the  requirements  for  the  doctoral  degree  at  the  University  of 
Illinois. 

*  Present  address:  Sandia  National  Laboratory,  Albuquerque, 
New  Mexico,  871S5.USA. 
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Figure  2  illustrates  that  current  interruption  is  a 
source  of  energy  in  a  plasma  focus.  The  lowest  trace 
is  the  main  bank  current  waveform  of  the  plasma 
focus.  The  sharp  reduction  in  the  current  by  a  factor 
of  two  at  pinch  time  indicates  [1]  that  one-half  the 
magnetic  energy  was  coupled  into  the  pinch  (i.e.  3.8  kJ). 
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spectra  of  ions  and  electrons  emitted  from  a  plasma 
focus  are  reported.  Also  the  scaling  of  the  primary 
electron  beam  current  and  energy  spectrum  with  the 
total  OPF  main  bank  current  at  pinch  time,  I  mb,  are 
reported  for  the  first  time  [4],  The  ions  and  electrons 
are  accelerated  in  the  directions  shown  in  Fig.  1.  The 
ion  beam  energy  spectra  were  measured  with  a  Faraday 


FIG. 4.  Schematic  view  of  electron  magnetic  spectrometer. 
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permit  direct  observation  of  fast  particles  emitted  from 
the  lower  end  of  the  pinch.  The  outer  diameter  of  the 
centre  electrode  is  S.08  cm  and  the  inner  diameter  of 
the  squirrel  cage  (outer  electrode)  is  10.2  cm.  The 
centre  electrode  extends  14  cm  beyond  a  lavite 
insulator  (26].  The  electrode  material  is  OFHC  copper; 
there  is  evidence  that  oxygen  released  from  regular 


0  7 


5S0  460  SCO 

OPF  MAIN  BANK  CURRENT  AT  PINCH  TIME  IN  KILOAMPS 
(LOG  SCALE) 


FIG. 5.  Plot  of  neutron  yield  versus  plasma  focus  current  at 
pinch  time.  Fillinf  pressure  was  S  tort  deuterium;  bank  roitape 
range  16-24  kK 
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used  previously  to  transport  electron  beams  (with 
v/y  <  1  and  vjy  *>  1)  with  less  than  20%  loss  of  beam 
charge  over  distances  of  ~  1  m  at  optimum .  ressure 
(31, 32],  The  optimum  pressure  is  defined  to  be  that 
pressure  at  which  the  electron  beam  charge  is  trans¬ 
ported  with  minimum  loss  (33,  34],  and  is  also  the 


return  current. 

The  electron  beam  Faraday  cup  that  was  built  is  a 
slightly  modified  [4]  Peilinen  [3]  Faraday  cup.  The 
cup  was  designed  to  measure  100-kA  currents,  has  an 
experimentally  determined  [4]  bandwidth  of  450  MHz, 
and  the  1/4  mil  mylar  plasma  current  filter  normally 
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3.  RESULTS 


Ed  0*V) 

FIG.6.  Ion  energy  spectra  from  four  shots  taken  on  25  June 
1980.  Charging  voltage  on  the  capacitor  bank  was  24  kV 
(11.5  kj);  pressure  of  deuterium  fill  gas  was  S.O  ton,  and 
neutron  yield  for  all  these  shots  was  within  10%  of  2.2  X109. 


The  results  of  the  energy  spectra  measurements  are 
shown  in  Figs  6—8.  In  Fig.6,  the  ion  spectral  data  for 
four  DPF  shots  are  shown  for  a  deuterium  fill  gas 
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their  plasma  focus  devices  is  optimal  (43). 

The  electron  magnetic  spectrometer  was  operated  as 
shown  in  Fig.4.  The  perpendicular  magnetic  field  was 
nominally  at  200  G  and  six  energy  channels  were  used, 
spanning  a  mean  energy  range  between  30  and  400  keV. 


FIG.  JO.  Typical  time-integrated  cltctron  energy  spectrum 
observed  with  the  electron  magnetic  spectrometer  /see  Fit-4) 
for  /mb"  370  kA.  Fillint  pressure  was  3  ton  deuterium;  bank 
voltage  was  17.5  kV  1/1/2)  CVx~b.l  kj). 
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of  the  electron  beam.  The  filter  had  to  be  replaced 
after  each  shot  and  so  a  gate  valve  was  inserted  into 
the  drift  tube  and  opened  after  the  plasma  focus  had 
been  conditioned. 
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signals  caused  by  these  electron  beams  which  have 
been  reported  in  the  past  (36-39], 

The  results  of  the  scaling  measurements  for  the 
primary  beam  current  are  shown  in  Fig.  13.  The  data 
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involving  L(V= q/C + LI)  estimated  for  If®  =  560  kA. 
To  compare  this  with  the  maximum  (Em)  energy  of  the 
particles  accelerated,  the  maximum-energy  channel 
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Waveform  a)  Is  delayed  40  tent  with  respect  to  oi,  aue  to 
Inherent  delay  of  photomultiplier  tube,  differences  In  cable 
length  and  llme-of-fltghi.  Filling  pressure  was  3  tort  deuterium; 
bank  voltage  was  23  kV. 
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if  a  few  simple  assumptions  are  made  about  the  target.  ilowtng-down  on  deuteron  beam  have  beer  Included.  77777 

If  the  target  is  assumed  to  be  the  cold  gas  ahead  of  the  '  Error  ban  are  at  in  iii)  and  thlt  letter  model  predim 

beam  then  it  is  independent  of  bank  parameters. 
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